Four-Step Evolution of Spin-Hall Conductance: Tight-Binding Electrons with Rashba 
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An intriguing magneto-transport property is demonstrated by tight-binding lattice electrons with 
Rashba spin-orbit coupling (SOC) in a magnetic field. With the flux strength <f> — 2n/N [N 
is an integer) and the Zeeman splitting fixed, when increasing the Rashba SOC A, the spin-Hall 
and charge-Hall conductances (SHC and CHC) undergo four-step evolutions: the SHC shows size- 
dependent resonances and jumps at three critical A c 's, and changes its sign at A c i and A C 3; while the 
CHC exhibits three quantum jumps by —Ne' 2 /h, +2Ne 2 /h and —Ne 2 /h. Such four-step evolutions 
are also reflected in topological characters and spin polarizations of edge states of a cylindrical 
system, and are robust against weak disorder. 

PACS numbers: 72.25.-b, 71.70.Ej, 73.43.Cd, 71.10.Ca 
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Introduction. — Recently, the spin-Hall effect (SHE), 
i.e., a generation of spin current perpendicular to an ap- 
plied electric field [H, [3, H, 0| , has shed new light on spin- 
tronics [B[ and provided novel techniques to manipulate 
spins in nanostructures. In contrast to the extrinsic SHE 
driven by spin-orbit (SO) impurity scattering [l], it is 
proposed that an intrinsic SHE exist in semiconductors 
with SO coupled bands @, Q ■ These proposals encour- 
aged the discovery of the SHE in GaAs semiconductor 
films and heterostructures and in metallic Al films 
and Pt strips [6]. In models with SO coupled bands, 
two-dimensional electron gas (2DEG) with Rashba spin- 
orbit coupling (SOC) Ml has the simplest form and is 
therefore most notable [a, lit Hoi. Illj ■ Meanwhile, tun- 
able Rashba SOC has been achieved via an external gate 
voltage on the top of asymmetric heterostructures [12j, 
and the Rashba SO field in quantum wells and semicon- 
ductors can also be measured optically [l3j |. 

In a clean 2DEG with parabolic dispersion and linear 
Rashba SOC, Sinova et al. predicted that the spin-Hall 
conductance (SHC) holds a universal value independent 
of SOC strength when both SO split bands are occu- 
pied [3(. It is now known that such an intrinsic SHC 
with only linear Rashba SOC might be destroyed by any 
amount of disorder [icj |. or be canceled completely by 
intrabranch contributions in the presence of a magnetic 
flux . In parallel, the SHC of 2DEG with linear Rashba 
SOC and Zeeman splitting in a magnetic field was cal- 
culated, and a resonant SHC was predicted when two 
Landau levels cross each other at the Fermi level 



11] 



In the presence of an underlying lattice potential, e.g. 
in metallic conductors like Al films and Pt strips Q , both 
parabolic dispersion and linear SOC should be modified 
and then incorporated into a lattice model which has 
been employed to study the effect of disorder on the 
SHE in the metallic regime Here, we begin to in- 
vestigate a lattice system of 2D tight-binding electrons 
(TBE) with Rashba SOC in a magnetic field. This model 



is also relevant to novel experimental systems such as 
untracold fermions in an optical lattice with an effec- 
tive SOC 14| and graphene with an intrinsic or Rashba 



soc [H Il6l |. We focus on magneto-transport properties, 
and have found that tuning the Rashba SOC strength 
generates novel four-step evolutions of the SHC and the 
charge-Hall conductance (CHC). Such bulk properties 
are also reflected in topological characters and spin po- 
larizations of edge states of a cylindrical system, and are 
robust against weak disorder. 

Formulation. — The model Hamiltonian of 2D TBE on 
a square lattice with Rashba SOC and a uniform perpen- 
dicular magnetic field B = (0,0, —B) is Q: 



H = -t^fe^cfo+H.c.] +A]T[z 



(ij) 



-ie 



o-yCi+3 + H.c. - h z y^ j (n lT - mi) (1) 



where cj 



(cj t , cj, ) are electron creation operators at 



site i, <j x and o~ y are Pauli matrices, the nearest-neighbor 
hopping integral t will be taken as the unit of energy, A 
is the Rashba SOC strength, and the Zeeman splitting 
parameter is h z = ^gfibB with g the Lande factor and 
fib the Bohr magneton. The magnetic flux per plaquette 
is <f) = J2n §n ~ ZTtBa 2 /<fio — 2ir/N with N an integer, 
a the lattice constant and </>o = hc/e the flux quantum. 
The Landau gauge A — (0, — Bx,0) and the correspond- 
ing periodical boundary conditions (PBCs) are adopted, 
and the magnetic unit cell has the size N X 1. 

After the numerical diagonalization of the Hamiltonian 
[Eq. (JTJ], the zero-temperature (T = 0) CHC is calcu- 



lated through the Kubo formula 



, „ x ie 2 h \ - \ - 
a CH (£) = — £ E 

e mk <Ee nk >E 

(mk|u 2 :|nk)(nk|wy|mk) — (mk|uy|nk)(nk|u x |mk) 



(2) 
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while the SHC at T = is given by 

Im{mk\J* s P iD \nk)(nk\v y \mk) 



eh \ - \ - 

<Ee nk >E 



(^mk ^nk) 



(3) 

where A = L x L is the area of this 2D system, E is 
the Fermi energy, is the corresponding eigenvalue 
of the eigenstate |mk) of mth Landau subband, and the 
summation over wave vector k is restricted to the mag- 
netic Brillouin zone (MBZ): -n/N < k x a < n/N and 
— 7r < k y a < 7r. The velocity operator is defined as 
v = i/h[H, R] (R is the position operator of electron) 
and the spin current operator as J x spm = h/A{v x , a z }. 
When E falling in energy gaps, we can rewrite ctch as 
(Tch(E) = e 2 /hJ2 Em< E Cm, where C m is the Chern num- 
ber 17[ of the mth totally filled Landau subband. 
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FIG. 1: (color online). The case with N = 4 and h = 0.4t. 

(a) The spin-Hall conductance ctsh versus the Rashba SOC 
parameter A for electron filling v = | and various lattice sizes. 

(b) The charge-Hall conductance itch versus A in the cases of 
(a), (c)-(j) The DOS for some A's in (a). The Chern numbers 
of subbands are also shown. 

An example with N = 4. — An overall picture of the 
CHC ctch and the SHC ct sh calculated by Eq. © and 
Eq. are shown in Fig. [T] with N — 4 (i.e., the flux 
strength cj> = j x 2ir), h z — OAt and various lattice sizes 
with L = 32—2048. We concentrate on the lowest Landau 
subbands and consider the electron filling v = i. 

In the case of A = [Fig. [U^c)] , the density of states 
(DOS) is symmetric about the Fermi energy E, and the 
lowest two Landau subbands (each totally-filled subband 
contributes | to v) are well separated, each carrying a 
Chern number +1. With A increasing from to l.Ot 
one sees a systematic four-step evolution of <tch and ctsh 
versus A; there are three critical A c 's at which both ctsh 
and itch exhibit jumps. 

When A increases from to A c i ~ OAAt, the lowest 
two Landau subbands approach each other, then merge 



together and form a pseudogap at A c i [Fig. [T^e)]; lt sh 
changes continuously from — le/8ir to larger negative val- 
ues [Fig. DJa)]; while ^ch = +le 2 /h nearly stays un- 
changed [Fig.[ljb)]. Here for a small lattice size (L = 32), 
<7sh and (tch both present divergence when A approaches 
A c i. With the lattice size increased (L — 64 — 512), the 
divergence is weakened accordingly; for L = 2048, osh 
approaches a finite value — 4.30e/87r at A c i, and ctch re- 
mains as +le 2 /h for < A < A c i- In the following, we 
focus on the data obtained with L — 2048. 

Increasing A across each A c , ctsh and ctch both exhibit 
sharp jumps: ctsh jumps from —4.30 to +4.87 (in units 
of e/87r) at A ci , from +8.71 to +6.32 at A c2 ~ 0.56i, 
and from +1.49 to —3.57 at A C 3 « 0.71t; ctch changes 
as +1 — » —3 — » +5 — * +1 (in units of e 2 /h). In inter- 
vals away from A c 's, ctsh varies continuously while ctch 
remains unchanged. The corresponding DOS [Fig. [lfc- 
j)] also points out that the lowest two Landau subbands 
approach, merge together and form a pseudogap at each 
A c , and then separate for three times. 

Mainly, such a four-step evolution of the SHC of TBE 
is distinct from the resonant SHC of 2DEG predicted by 
Shen et al. [ll| in four aspects: in 2DEG, two Landau 
levels cross each other at the Fermi level only once and 
produce one A c , while for TBE the two Landau subbands 
touch successively three times and results in three A c 's; 
at a A c , the SHC of 2DEG diverges at T = 0, while the 
SHC of TBE converges to finite values in the thermody- 
namic limit (L -> oo) at T = 0; the SHC of 2DEG does 
not change its sign while the SHC of TBE changes its 
sign at Ad and A C 3; furthermore, the CHC of 2DEG is 
unaffected when tuning A, but the CHC of TBE presents 
three successive quantum jumps. 
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FIG. 2: (color online), ush versus the Rashba SOC parameter 
A in various cases, ctch (in units of e 2 /h) of each evolution 
step is also shown. 

Cases with weaker magnetic fields. — The above four- 
step evolutions have also been verified by further numer- 
ical calculations of the cases with TV = 4 — 16, h z = 
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0.05i — OAt, and various v's (with odd number of totally 
filled Landau subbands), as illustrated by four examples 
in Fig. EU For N = 4, h z = 0.2t and v = 1/8 [Fig. C^a)], 
ctsh shows behaviors similar to that in Fig. [TJa), while 
with smaller A c 's and narrower transition regions (i.e., 
smaller A c3 - A cl ); for N = 6 and N = 8 [Fig. db-d)], 
the transition regions are narrower than the case with 
N = 4. Meanwhile, the quantized CHC also exhibits 
three jumps by -Ne 2 /h, +2Ne 2 /h and -Ne 2 /h. 

In brief, the larger iV's, the significantly narrower are 
the transition regions (A c i < A < A C 3). However, the 
positive values in the transition regions are much larger, 
and the total weights of positive part of osh (he. the 
integral from A c i to A C 3) possessing the same order of 
magnitude, are respectively 0.96, 1.69, 0.99 and 1.12 in 
the four cases of Fig. [2] [Note that the weight is 1.22 for 
the case in Fig.QJa).] 




FIG. 3: (color online). Lowest two subbands and intermediate 
edge states [shown as thick (red) lines] of a cylinder of the size 
128 x oo (OBC in x direction and PBC in y direction) with 
N — 4, h z — 0.4 and various A's. 

Edge states in a cylindrical system with N — 4. — An 
alternative way to reveal the distinctions among four evo- 
lution steps is to calculate the edge states of the system 
on a cylinder. These edge states reflect the topological 
character of the corresponding bulk state Tj^, 201 ]. Just 
recently, spin-filtered edge states have been considered 
for a graphene cylinder with an intrinsic SOC [H[ (in 
a two-component Haldane model [2l| ) or Zeeman split- 
ting , a quantum SHE arising from helical edge states 
have been proposed and experimentally verified soon in 
HgTe quantum wells [22j , and edge states have also been 
employed to characterize topological band insulators and 
chiral spin liquids [23} . Now as an illustration, we take a 
cylinder of square lattice of the size 128 x 00 and apply 



open boundary condition (OBC) in x direction and PBC 
in y direction. 

Chern numbers of bulk Landau subbands are inti- 
mately related to the winding numbers of the correspond- 
ing edge states 20]. For A ci < A < X c2 , there is one edge 
state winding three times from the upper subband to the 
lower one then back to the upper one [a thick (red) line 
in Fig. [3jb)] which corresponds to a Chern number —3 of 
the lower subband. For X c2 < A < A C 3, there is one edge 
state winding five times from the lower subband to the 
upper one then back to the lower one [Fig. [3jd)] which 
corresponds to a Chern number +5 of the lower subband. 
While for < A < A c i or A C 3 < A < l.Ot (not shown in 
Fig. [3]), there is another edge state winding only once 
from the lower subband to the upper one then back to 
the lower one which corresponds to a Chern number +1 
of the lower subband. 

The continuum spectrum of this cylinder also gives fur- 
ther descriptions about the jumps of the bulk CHC. In- 
creasing A across A c i [Fig. [3ja)] or A C 3 [Fig. [3je)], the 
lowest two subbands touch at four points simultaneously 
in fc-space and a Chern number —4 is transfered from 
the upper subband to the lower one; while across A C 2 
[Fig.(3Jc)], the lowest two subbands touch at eight points 
simultaneously in fc-space and a Chern number +8 is 
transferred between them. Such a correspondence be- 
tween transferred chern numbers and touching points in 
/c-space has also been verified for iV = 5 — 8. 



0.08 

z 06 
2 

^0.02 
J 00 



(c) 




(g) 



0.06 

I 03 

■c 0.00 
o 

Q--0.03 
0.04 
§ .03 

H 02 
0.01 
£ 0.00 
0.1 



(b) 



(d) 



1 < f ) 



(h) 



Lattice Site Index 



Lattice Site Index 



FIG. 4: (color online). Spin polarization P z (in units of h/2) 
versus the lattice site index in x direction for the edge states 
of the cylindrical system in Fig. [3] 



In addition, the spin polarization carried by the 
edge states can be computed explicitly as P^ k (i) = 
ft/2(mk c\ a z Ci\m\L) with i the lattice site index in x di- 
rection 16|. In Fig. 21 we plot the spin polarization P z 
of some edge states of the above cylindrical system. If A 
takes a value far away from A c 's, P z takes prominently 
large values near the left or the right edge and is almost 
zero in the intermediate region [Fig. Ufa), (c) and (f-h)]; 
but if A takes a value close to A c 's, P z fluctuates strongly 
between two edges [Fig. |3Jb), (d) and (e)]. Note that for 
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a fixed k y , the dominantly positive peak of P z moves to 
another edge when A varies from OAt to 0.5t. And for 
edge states of A = 0.8t [Fig.QJh)], P z takes prominently 
negative values near edges. 

Presence of disorder. — We add a term iDjcjci Q 
into the Hamiltonian [Eq. (fl])] to account for the effects 
of nonmagnetic disorder, Wi being a random potential 
uniformly distributed between [— W/2,W/2]. 



and graphene with an intrinsic or Rashba SOC [SGI. 
And spin polarizations of edge states should be observ- 
able in a four-terminal experimental setup 15, ljj ]. 




FIG. 5: (color online). <tsh versus A in the case with N = 4, 
h z — OAt and various disorder strength W's (100 random- 
potential configurations of the size 8x8). The inset shows 
the evolution of <tsh versus W at A = 0.5t. 

For N = 4, the adopted 100 random-potential config- 
urations are of the size 8x8 (such a super unit cell is 
commensurate with the magnetic unit cell in the absence 
of disorder), and the total lattice is of the size 32 x 32. 
It can be seen that weak disorder (W < 0.5i) does not 
smear out the overall four-step evolution of the SHC. For 
stronger disorder (W = 2M), the SHC does not show 
resonance anymore near A c i or A C 3, and takes positive 
values in an enlarged interval while the peak is dimin- 
ished into a hump. 

Summary and discussion. — An appealing evolution of 
magneto-transport property has been demonstrated by 
TBE with Rashba SOC in a magnetic field: (i) with the 
flux strength <fi = 2ir/N and the Zeeman splitting fixed, 
when increasing the Rashba SOC A from 0, four-step evo- 
lutions of the SHC and CHC have been observed; (ii) at 
three A c 's, the SHC shows size-dependent resonances and 
jumps, and changes its sign at A c i and A C 3; (iii) mean- 
while, the quantized CHC shows three successive jumps 
by -Ne 2 /h, +2Ne 2 /h and -Ne 2 /h; (iv) for smaller 0's, 
the total weights of positive part of SHC have the same 
order of magnitude although the transition regions are 
significantly narrower; (v) edge states of a cylindrical 
system reflect such bulk properties; (vi) this four-step 
evolution is robust against weak disorder. 

Such a four-step evolution of SHC is expected to occur 
in 2D electron systems with a lattice potential, a mecha- 
nism of SOC or SO scattering, and an external magnetic 
field. Some candidate experimental systems are: metallic 
conductors such as Al films and Pt strips Q, untracold 
fcrmions in an optical lattice with an effective SOC 14}, 
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